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We investigate the long-time behavior of a majority rule opinion dynamics model in finite spatial 
dimensions. Each site of the system is endowed with a two-state spin variable that evolves by 
majority rule. In a single update event, a group of spins with a fixed (odd) size is specified and 
all members of the group adopt the local majority state. Repeated application of this update 
step leads to a coarsening mosaic of spin domains and ultimate consensus in a finite system. The 
approach to consensus is governed by two disparate time scales, with the longer time scale arising 
from realizations in which spins organize into coherent single-opinion bands. The consequences of 
this geometrical organization on the long-time kinetics are explored. 
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I. INTRODUCTION 

The majority rule model (MR) is a simple descrip- 
tion for consensus formation in an interacting population. 
The model consists of spins (opinions) that are fixed 
on lattice sites, and each spin can assume the states -1-1 or 
— 1, corresponding to two opposite opinions. Spins evolve 
by the following two steps: first, pick a group of spins of 
fixed odd size G; second, all the spins in this group adopt 
the state of the local group majority (Fig. These two 
steps are repeated until a final consensus is necessarily 
reached. Our goal is to understand basic properties of 
this approach to consensus in finite spatial dimensions. 





FIG. 1; Illustration of a single majority rule update step for 
the 5-site von Neumann neighborhood on the square lattice. 

A general form of this majority rule dynamics was 
introduced by Galam 0| in which a variable number 
of groups of arbitrary size are formed simultaneously 
and then majority rule is simultaneously applied to each 
group. Our implementation of majority rule, in which 
only a single small group is updated at each time step, 
allows for considerable analytical progress in the mean- 
field limit 0, 1^ and also makes it convenient to simulate 
the model, especially in high dimensions. 

In a previous study of the MR model , it was shown 
that the average time until consensus is reached is pro- 
portional to the logarithm of the number of spins N in the 
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system in the mean-field limit. On the other hand, for 
finite dimensions, numerical simulations suggested that 
the most probable consensus time grows as a power law 
in N , with an exponent that decreases as the spatial di- 
mension increases. Mean-field behavior was not repro- 
duced even in four dimensions, indicating a still larger 
value for the upper critical dimension of the MR model. 

In this article, we focus on the MR model in finite spa- 
tial dimensions. The questions that we will investigate 
are: What is the geometry of single-opinion domains? 
How long does it take to reach consensus? How do basic 
system parameters affect the consensus time? We find 
that the probability distribution for the consensus time 
involves two very different time scales when the spatial di- 
mension is greater than one. The longer time scale arises 
from configurations in which opposite-opinion domains 
organize into coherent geometries - stripes in two dimen- 
sions, slabs in three dimension, etc. While the probability 
for the system to reach such a coherent state decreases as 
the spatial dimension is increased - approximately 33% in 
two dimensions and 8% in three dimensions - we believe 
that this probability remains non-zero in all finite spatial 
dimension. More importantly, the time needed to reach 
final consensus from these coherent states is extremely 
long. These configurations therefore give the dominant 
contribution to the mean consensus time. 

To put our results in context, it is instructive to com- 
pare the MR model with two fundamental kinetic spin 
models, namely, the voter model (VM) j^, and the ki- 
netic Ising model with zero-temperature Glauber kinet- 
ics (IG) 0. The VM describes consensus formation in a 
population of individuals with zero self confidence. In an 
update step of the VM, a spin is selected at random and 
it blindly adopts the state of a randomly-selected neigh- 
bor. This step is repeated until consensus is necessarily 
reached. Because of the underlying linearity of the VM 
spin-flip rate on the number of anti-aligned nearest neigh- 
bors, the VM is exactly soluble in all spatial dimensions 
^1^0. In particular, for an A^-spin system in d dimen- 
sions with zero initial magnetization, the consensus time 
scales as N for d > 2, as A^^lniV in d = 2 (the critical 
dimension of the VM), and as N'^ in d— 1. Because the 
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average magnetization is conserved, the probability that 
the system eventually ends with all + spins equals the 
initial density of + spins in all spatial dimensions. 

In contrast, the zero-temperature kinetics of the IG 
model obeys a form of majority rule. In the update step, 
a flippable spin (those with zero or positive energy) is 
picked at random and it adopts the state of the majority 
in its interaction neighborhood. In the case of a tie in 
the neighborhood state (which can happen on bipartite 
lattices), the selected spin flips with probability 1/2. This 
elemental update step is repeated until no flippable spins 
are left. At early times, coarsening domains form whose 
typical length scale grows as i^/^ due to an underlying 
diffusive dynamics Q. 

The primary operational difference between the IG and 
MR models is that in the latter all the spins within the 
neighborhood flip, a feature that also occurs in Galam's 
model PI and also in the Sznajd model 9] of social in- 
fluence, where a small group that is in consensus can 
influence other spins at the periphery of the group. This 
distinction in the update rule has fundamental conse- 
quences. In the IG model, infinitely long-lived metastable 
states can occur that consist of perfectly flat interfaces 
in cZ = 2, or states where all interfaces have zero net cur- 
vature for d > 3 Ref. contrast, consensus is the 
only possible final state in the MR model. Nevertheless, 
both the zero-temperature IG model and the MR model 
have anomalous kinetics because of the existence of very 
long-lived transient states. 

In Sec. II, we present simulation results for the anoma- 
lous behavior of the consensus time distribution and the 
two basic controlling time scales. Then in Sec. Ill, we dis- 
cuss the role of the long-lived coherent states that dom- 
inate the asymptotic tail of the consensus time distri- 
bution. A qualitative argument for the lifetime of these 
states is given in Sec. IV. We conclude in Sec. V. 



II. CONSENSUS TIME DISTRIBUTION 

We first simulate the distribution of times until con- 
sensus is reached on finite-dimensional hypercubic lat- 
tices with periodic boundary conditions. Typically, we 
initialize each realization of the system to contain equal 
numbers of + and — spins. We choose the group size to be 
G = 3 and construct the group by selecting a spin at ran- 
dom and then randomly picking two out of its 2d nearest 
neighbors. This definition for a group has the advantages 
of computational simplicity and a dimension-independent 
group size. Other definitions for a group, such as the von 
Neumann neighborhood of Fig. (the initial site plus its 
2d nearest neighbors; group size G = 2d + 1), lead to 
qualitatively similar results. 

We then evolve each realization according to MR ki- 
netics until consensus is reached. The quantities that 
we focus on are: (i) the distribution of consensus times, 
P/v(t), in an TV-spin system with zero initial magnetiza- 
tion, and (ii) the probability for a realization to reach a 



stripe or a slab state, S'Ar(m) (to be defined below), as a 
function of N and the initial magnetization m. 
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FIG. 2: Distribution of consensus times, Pjv(t), versus t (in 
units of Monte Carlo steps per spin) for an A''-spin system 
starting from a random zero-magnetization state. Shown are 
results for (top to bottom): a 50^ square lattice (A^ — 2500), 
a 14'^ cubic lattice (A'^ = 2744), and a four- dimensional 
hypercubic lattice (A'^ — 2401). Data are all based on 10^ 
realizations. 
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The consensus time distributions PN{t) for spatial di- 
mensions 2, 3, and 4 are shown in Fig. |21 It is evident 
that in two and three dimensions, PN{t) is character- 
ized by two time scales - the most probable consensus 
time, corresponding to the peak of the distribution, and 
a much longer time scale associated with the asymptotic 
exponential decay. In four dimensions, there is a change 
in the slope of the asymptotic tail of PN{t) for t > 400, 
suggesting the possibility that the asymptotic kinetics 
involves yet a third time scale. 
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FIG. 3: Double logarithmic plot of Mk{N), the fc"' reduced 
moments of the consensus time distribution, versus the total 
number of spins iV for A; = 1 (o), it = 2 (□), fc = 4 (A), fc = 8 
in d = 2 (top) and d — 3 (bottom). Data are based on 
10^ realizations in d = 2 and 4 x 10^ realizations in d = 3. 
The lines are least-squares linear fits to the large- data. 



we now apply two methods to estimate this longer time 
scale. First, we consider the reduced moments of the 
consensus time distribution 



MkiN) = ((t(iV))'=)i/^- 



t'' PN{t)dt 



l/k 



(1) 



As suggested by the data in Fig. 13 if the long-time tail 
of consensus time distribution has a simple exponential 
decay of the form e~*/'^*^^^ at long times, then all the 
reduced moments would asymptotically scale as t{N), 
with subdominant corrections that become smaller as k 
increases. This trend is illustrated in Fig. where Mk 
is plotted as a function of TV for various values of k. In 
d = 2, each grows as a power law in TV for large k, 
but with a slightly different apparent exponent. Least- 
squares fits to the data give the following exponents in 
d = 2: 1.64 for fc = 1 (o), 1.73 for fc = 2 (□), 1.75 for 
k = A (A), and 1.75 for /c = 8 (+). From this hmiting 
large-fc value of this exponent we can then infer the N 
dependence of r. 

For d = 3, the behavior is qualitatively similar, except 
that there is a large disparity in the exponents for 
ioi k — 1 and for A: > 1. Linear fits to the data now 
give the exponent values 1.18 for fc = 1 (o), 1.64 for 
fc = 2 (□), 1.77 for fc = 4 (A), and 1.78 for fc = 8 {+). 
However, there is a perceptible downward curvature in 
the dependence of Mk on N for large fc, so that linear 
fits are inadequate to determine the A'^-dependence of r 
accurately. 
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FIG. 4: Double logarithmic plot of t{N) (in units of Monte 
Carlo steps per spin) versus for d = 2 (o) and d = 3 (A). 
The lines are the best fits to the last few data points. 



From these data, we find that the most probable con- 
sensus time, ^mp, scales with TV as imp ~ TV", with 
a ~ 1.24, 0.72, and 0.56 for spatial dimensions 2, 3, and 
4, respectively. These values are identical to those ob- 
tained previously in Ref . ; these were based on smaller- 
scale simulations in which realizations where the consen- 
sus time exceeded a (large) preset limit were terminated. 

On the other hand, the asymptotic decay of P/v(t) is 
clearly governed by a much longer characteristic time and 



Our second analysis method is simply to measure the 
slope of the exponential tail of Pm {t) directly for differ- 
ent values of TV and thereby determine t(TV). To do this, 
we first make a first estimate for t by finding the slope 
in the region that is visually most linear. Then we refine 
this estimate by computing the slopes in the systematic 
ranges (r/2,2r), (T/2,3r), (T/2,4r), etc., and using the 
range where a linear fit has the highest correlation co- 
efficient. In the resulting data for t versus TV, there is 
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now small and systematic downward curvature (Fig.^. 
By dropping the first four data points one-by-one and 
then performing linear fits to the remaining data, the 
local slope decreases from 1.746 to 1.719 in two dimen- 
sions. In three dimensions, there is a larger decrease in 
the local slope from 1.832 to 1.709 as the first six points 
are deleted. Extrapolating this local slope to oo, 
we obtain the estimates i/ = 1.7 ± 0.04 in d = 2 and 
1/ = 1.5 ± 0.1 in d = 3 in the relation r ~ A^". The er- 
ror bars are a subjective guess of the uncertainty in the 
extrapolation. 



III. ANOMALOUS COARSENING AND 
LONG-LIVED COHERENT STATES 

The main result of the above analysis is that the aver- 
age consensus time is much larger than the dependence of 
N^^'^ that would arise if domain coarsening were entirely 
governed by diffusive dynamics. By observing the evolu- 
tion of many realizations of the system, it is clear that the 
asymptotic tail of the consensus time distribution arises 
from situations where the -|- and — spins organize into 
spatially coherent and long-lived states that consist of 
relatively flat stripes in two dimensions (Fig.[SJ), slabs in 
three dimensions (Fig. |HJl, and analogously (we believe) 
in higher dimensions. The existence of these states is one 
of the most surprising feature of the MR model. In spite 
of the isotropy of the MR interaction, the long-lived tran- 
sient states arise and spontaneously break this symmetry. 
Once the system reaches such a state, further evolution 
proceeds extremely slowly, as we shall discuss below. 

To develop intuition for these coherent states, we show 
in Fig. |S1 a set of snapshots of a 50 x 50 system that 
happens to evolve to a stripe. After a few time steps, 
the lattice-scale granularity of the random initial state 
has disappeared due to the effective surface tension in 
the majority rule dynamics. After this early-time tran- 
sient, the subsequent evolution qualitatively resembles 
the coarsening of a spin system with non-conserved order- 
parameter kinetics. However, domains tend to develop a 
stringy morphology, a feature that promotes the forma- 
tion of stripes that span the system. For the realiza- 
tion shown, a clearly resolved stripe emerges by 100 time 
steps, while ultimate consensus is achieved when 1850 
time steps have elapsed; notice that a time of 1850 steps 
is relatively early in the asymptotic tail of P/v(i) in Fig.|2] 

In spite of the anomalous long-time kinetics of the MR 
model, the early-time coarsening is diffusive in nature. 
To determine the growth of the typical domain length 
scale at early times, we studied the time evolution of the 
two-spin correlation function. We took this correlation 
function at different times and found the length rescaling 
that gave the best data collapse. We thus found that 
the appropriate rescaling the correlation function is by 
a length scale that is proportional to t^^"^. We therefore 
conclude that the early-time coarsening in the MR model 
is characterized by a length scale that grows as t-^^"^. 




FIG. 5: Snapshots of a 50^ system at t = 0, 1, 5, 20, 80, & 200. 

A phenomenon analogous to stripe formation occurs in 
three dimensions, where long-lived states arise that con- 
sist of two relatively flat slabs of oppositely-oriented spins 
(Fig.|BJ. For the example shown from a 20^ lattice, a slab 
state forms around 150 time steps, while final consensus 
does not occur until 3800 time steps have elapsed. 



X 




FIG. 6: One of the two interfaces of the slab state on a cubic 
lattice of linear dimension 20. Coordinates are in units of 
lattice spacing. 

To verify that stripe states actually govern the asymp- 
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totic tail of PN{t) two dimensions, we also study the evo- 
lution of a synthetic system with an ordered initial state 
that consists of two straight stripes, with half the spins 
+ and half the spins — . The long-time tail of the con- 
sensus time distribution for this special initial condition 
follows a single exponentially decaying function, as shown 
in Fig. [7| Also shown in this figure is the corresponding 
distribution for a system of the same size with a random 
zero-magnetization initial condition. The coincidence of 
the slopes in the tails of these two distributions shows 
that stripe states control the long-time evolution of ran- 
dom zero-magnetization initial condition systems. 

Because of the crucial role that spatially coherent 
states play in the MR model, we also study the prob- 
ability SNim) that a randomly-prepared A'^-spin system 
in d dimensions with initial magnetization m evolves to 
such a state. We use two independent methods to mea- 
sure SN{m). One is based on simply counting the frac- 
tion of realizations whose consensus time lies within the 
asymptotic tail of the consensus time distribution. For 
example, for the data from the 50 x 50 system in Fig. |21 
the tail region corresponds to a consensus time t greater 
than 600. Thus for this system size all realizations with 
t > 600 are counted as reaching a coherent state. 




5000 10000 15000 20000 
t 

FIG. 7: Consensus time distribution for a 50^ square lattice 
for a two-stripe (dashed) and random initial condition (solid). 
All data are based on 10^ realizations. Time is in units of 
Monte Carlo steps per spin. 

Alternatively, we investigate correlation functions that 
are engineered to detect stripe states. For d = 2, we 
consider the following correlation functions for two spins 
that are located a distance L/2 apart: 

C,{t) = ^{s{x,y,t)s{x±L/2,y,t)), 
Cy{t) = ]^{s{x,y,t)s{x,y±L/2,t)). 

For both a random state and for consensus, these corre- 
lation functions equal zero. Conversely, for an ordered 



two-stripe state with stripes of width L/2 parallel to the 
X-axis, Cx = -1-1 and Cy = —1, and vice versa for stripes 
parallel to the y-axis. We therefore posit that a stripe 
state arises if the correlation function in the direction(s) 
parallel to the stripe is greater than a threshold value, 
while the correlation function perpendicular to the stripe 
is less than the negative of this threshold value. We arbi- 
trarily choose the threshold to equal 0.5, but our results 
for large N depend only weakly (S'Ar(m = 0) varies by 
< 10%) on the threshold value when it is in the range 0.3 
- 0.7. The results given below are based on the threshold 
set to 0.5. 

0.4 1 ' , ' , ' , ' , ' 1 
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FIG. 8; Probability to reach a stripe state versus number of 
spins A for d = 2 (o) and d = 3 (□) using the threshold value 
0.5 (see text). Data are based on 10^ realizations. 



We find that the stripe/slab probability Sj^(m — 0) 
grows quickly for small A^ and then saturates to a non- 
zero value that is close to 0.33 in d = 2 and 0.08 for d = 3 
(Fig. IHl) . The stripe probability in two dimensions is very 
close to that found previously in the zero-temperature 
evolution of the Ising model with Glauber kinetics [Toj . 
Note also that as the initial magnetization m is moved 
away zero, S^ijn) quickly decays to zero (Fig. |5J). This 
simply reflects the fact that if one phase is initially below 
the percolation threshold, there is a very small possibility 
for minority phase droplets to merge and form a stripe 
that spans the system. 

We can qualitatively understand the dimension depen- 
dence of the probability to reach a stripe state by the 
following rough argument (see also Rcf. T^). For sim- 
plicity, we first discuss the case of two dimensions with 
the random zero-magnetization initial condition. Con- 
sider a large system of linear dimension L and cut it into 
four equal subsquares of linear size L/2. The final state 
in each of these subsquares is reached more quickly than 
that of the entire system. We now make the plausible 
assumption, based on observations of many realizations 
of the system, that each subsquare independently reaches 
consensus. Then out of the 2'* possible configurations of 
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FIG. 10: Final magnetization as a function of initial mag- 
netization m, on the square lattice for linear dimension L = 10 
FIG. 9: Probability to reach a stripe state versus initial mag- (^^^^^^^ ^^^^^^^ 20 (dashed), 30 (dot-dash), and 50 (solid), 
netization m for a 50 system. Data are based on 10"" realiza- 
tions. 



these subsquares, only the following arrangements 




where the -I- and — symbols refer to the final state of 
each subsquare, correspond to a stripe state of the L x L 
system. This argument then suggests that SN{m = 0) = 
4/2^ = 1/4. 

This coarse-graining argument straightforwardly gen- 
eralizes to higher spatial dimensions. On the cubic lat- 
tice, we divide a,n L x L x L cube into eight subcubes of 
linear dimension L/2. If these subcubes each indepen- 
dently reach consensus, then a slab state on the original 
cube (consisting of two slabs of oppositely oriented spins, 
each of size L x L x L/2) can be achieved in six possible 
ways. The probability of reaching a slab state is therefore 
6/2^ Ri 0.047. In d dimensions, this same line of reason- 
ing gives S'Ar(m = 0) = 2d/2^ . While our argument is 
crude, the resulting numerical values for Sn qualitatively 
mirror the corresponding estimates from simulations. 

Our approach also helps explain why stripe states 
quickly disappear when the initial magnetization is non- 
zero. As an example, for initial magnetization 0.08, we 
find by numerical simulations that the probability that a 
25 X 25 system eventually ends with all spins -|- is 0.88. 
Now employing the above coarse-graining argument for a 
50 X 50 system, the four 25 x 25 subsquares will each reach 
-|- consensus with probability 0.88 and — consensus with 
probability 0.12. Then the probabihty for the 50 x 50 
system to reach a stripe state is Sn ~ 4(0.88)^ (0.12)^ sa 
0.0446. This is very close (probably fortuitously) with 
our numerical result of SN{m = 0.08) ~ 0.0498. 

Another important aspect of the evolution to the final 
state is the dependence of the final magnetization on the 
initial magnetization. Since the system always reaches 
consensus, the final magnetization is simply the dif- 



ference in the probabilities that the systems ends with all 
spins -|- and all spins — . On the square lattice, we find 
that the curve of the final magnetization versus the ini- 
tial magnetization approaches a step function as — + cxd 
(Fig. I10|l . Thus any initial bias pre-determines the final 
state of the system in the thermodynamic limit. This 
step-function behavior is in contrast to the behavior in 
one dimension, where the final magnetization curve re- 
mains non-singular as A^ ^ cxd 0. Finally, it is worth 
noting that nif equals the initial magnetization m for 
the voter model in all spatial dimensions j^]; there is no 
tyranny of the majority in the voter model. 



IV. LIFETIME OF STRIPE STATES 

Once a stripe state is formed, the evolution to ultimate 
consensus is controlled by the time required for the two 
interfaces that define the stripe to meet and annihilate. 
In one dimension, it is easy to see that each isolated in- 
terface between -f- and — spins moves by free diffusion. 
When two interfaces approach to nearest-neighbor sep- 
aration they necessarily annihilate. (Note that in the 
kinetic Ising model, two nearest-neighbor interfaces can 
annihilate, with probability 1/2, or recede by one lat- 
tice spacing, also with probability 1/2, in a single update 
step.) Therefore the time for the last two domain walls 
to annihilate is proportional to L^, where L is the linear 
dimension of the system. Further, because the system is 
controlled by the meeting of two random walks on a finite 
ring, the consensus time distribution has an exponential 
decay of the form e~*^^ 0. 

In two and three dimensions, the interfaces between 
stripes are quite smooth (Figs. |S1 & and the scahng 
of the interface width on the transverse dimension of the 
system appears to be in the Edwards- Wilkinson univer- 
sality class !Tl|. We verified the smoothness of the in- 
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terface by preparing a system of linear size £ x 
in d dimensions, with £ 3> i, in which all spins in the 
region [0, £/2] are initially in the — state, and all spins 
in the region [£/2,>C] are initially in the + state. In two 
dimensions, the width w of the interface initially grows 
slowly in time and eventually saturates to a value that 
approximately scales 'as w L^/^. In three dimensions, 
the growth of the width is even slower and the satura- 
tion value of the width is consistent with a logarithmic 
dependence on L. 

Thus it is the diffusion of the interface as a whole rather 
than fluctuations in the interface shape that determines 
the lifetime of the stripe state. Since the interfaces are 
typically separated by a distance of order L when they 
are first formed, the lifetime T of the stripe state should 
therefore given by 



T . 



D{Ly 



(2) 



where D{L) is the diffusion coefficient of a single interface 
with transverse dimension L. 

We may obtain a simple albeit rough estimate for this 
diffusion coefficient by treating each site on the inter- 
face as an independent random walk [lol . For a d- 
dimensional system, a smooth interface contains of the 
order of L''"^ sites. In a single time step, each inter- 
face site will randomly move by ±1 perpendicular to 
the interface. Hence if each site is independent, the 
center of mass of the interface will move by a distance 
vZ/^^/L'*^^ ~ ^-(^-1)72 in one time step. As a result, 
the diffusion coefficient of the interface D{L) scales as 

We tested this prediction by simulation by following 
the evolution of a single interface in a long strip (or slab) 
geometry with transverse dimension L in which all the 
spins on the right half are set to -1-1 and all the spins 
on the left half are set to —1. We then let the spins 
evolve by majority rule dynamics. After a short transient 
that lasts of the order of one time step, we observe that 
the interface moves diffusively, with a diffusion coefficient 
that scales approximately as in two dimensions and 
as in three dimensions. Given the crudeness of the 
above random walk argument, it is surprising that the 
simulation results agree quite well with the prediction 
D{L) - L-^'^-^\ 

From this scaling of the diffusion coefficient on the 
transverse linear dimension L, Eq. (0) then gives a con- 
sensus time T that scales as T ~ L'^^^. Equivalently, in 
terms of the total number of spins N — L"^, the depen- 
dence is r ^ ]\[id+i)/d^ However, this prediction is only 
qualitatively consistent with the exponent values of 1.7 
for d = 2 and 1.5 for d = 3 that were obtained from direct 
numerical simulations of the consensus time distribution. 
We do not have an explanation for this discrepancy. 



V. SUMMARY AND DISCUSSION 

We studied the time evolution of the majority rule 
(MR) model for finite-dimensional systems. One of our 
main results is that the approach to consensus in an ini- 
tially unbiased system is surprisingly complex. Before 
ultimate consensus is reached, a non-trivial fraction of 
all realizations falls into coherent metastable states that 
consist of stripes in two dimensions and slabs in three di- 
mensions. We anticipate that analogous coherent states 
arise in higher spatial dimensions. The interfaces be- 
tween domains in these coherent states are quite smooth 
and reflect the strong surface tension in the majority rule 
dynamics. 

Due to these coherent states, the time to reach consen- 
sus is anomalously long and is controlled by a diffusion 
process that brings two interfaces close enough that they 
can annihilate. The characteristic time scale for this an- 
nihilation is much longer than the most probable time to 
reach consensus. The fraction of realizations that reach 
these long-lived states decreases as d increases, but their 
role appears to be dominant in the asymptotic kinet- 
ics. When the initial magnetization is non-zero, however, 
these long-lived states quickly disappear. As a result, the 
consensus time distribution has a single peak and there is 
no long-time tail. Furthermore, the time until consensus 
grows only logarithmically in the system size fFig. Ill|l . 
Thus an initial bias in the density of spins is a decisive 
influence in the long-time behavior of the system. 
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FIG. 11: Consensus time (in units of Monte Carlo steps per 
spin) versus A'^ in d = 2 for initial magnetization m = 0.2. 

We gave a crude coarse-graining argument to estimate 
the probability to reach a coherent state as a function 
of the spatial dimension. This approach qualitatively ex- 
plained the behavior of the probability to reach such a 
state as a function of the spatial dimension d and the 
initial magnetization. 

Finally, we suggest several directions for further study. 
First, it would be worthwhile to determine the value of 
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the upper critical dimension of the MR model. An exact 
analysis of this model on the complete graph, where all 
spins are nearest neighbors of each other, showed that 
the mean consensus time grows as In TV Q ■ On the other 
hand, the simulations presented here and in suggested 
that the mean consensus time grows as a power law in N 
for spatial dimensions 1, 2, 3, and 4. These two facts sug- 
gest that the upper critical dimension of the MR model 
is greater than 4. It would be worthwhile to have a the- 
oretical understanding for the apparently large value of 
the upper critical dimension. 

We also believe it will be fruitful to study simple ex- 
tensions of the MR model with more stringent conditions 
for achieving consensus. One example is to have a higher 
threshold than simple majority before the opinion of a 
group is swayed. While a higher threshold will obviously 
slow the dynamics, it should be interesting to investigate 
whether this modification leads to different scaling prop- 
erties for the mean consensus time and the distribution 
of consensus times. 

A more intriguing generalization arises when each spin 
has more than two opinions, where we anticipate new 
types of dynamical behavior. With more than two states, 
the possibility of a dynamically-stable steady state that 
consists of coalescing and coexisting multiple opinion 
groups was discussed in the framework of the "stochastic 
seceder" model |13J. In the context of our majority rule 
model, there obviously will be slower dynamics because 
it may be possible to have a group with no local majority. 



but only a local plurality. Such a group would not evolve 
according to the majority rule dynamics. Thus configu- 
rations in which there is no majority in each group rep- 
resent another absorbing state for the dynamics. In the 
mean- field limit, we find that such a system never reaches 
this frustrated state, as the corresponding fixed^oint of 
equal concentrations of all species is unstable [l^]. In- 
stead, for a system with more than two opinion states, the 
time to reach ultimate consensus is merely increased by 
a multiplicative factor compared to the two-opinion MR 
model. However, for finite spatial dimensions, the exis- 
tence of more than two opinions appears to have a more 
significant effect on the long-time behavior that depends 
fundamentally on the interplay between the group size 
and the number of states. When there are many distinct 
local majorities in the initial state the group dynamics 
has a primarily diffusive character. However, when there 
is of the order of one local majority, the opinion of this 
group quickly overtakes the entire system. 
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